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ABSTRACT 



The system studied is a surface-to-air guided missile whose 
aerodynamic and control characteristics can be represented by a set 
of linear differential equations. The objects of the study are to select 
a set of control system gains which will give specified performance of 
the missile within the roll and yaw subsystems and at the same time 
to minimize the coupling effects whereby motion about one of the axes 
causes dynamic response about the other axis. 

The system is studied analytically using the mathematics of 
control system synthesis and design. Analytic investigations are 
paralleled by system simulation on a REAC analogue computer. 

Using the above methods, gains are selected which give 
specified response in the roll stabilization and yaw control subsystems 
and suitably damped response due to coupled motions between systems. 
The study also investigates the possibilities of compensation to 
attenuate or eliminate coupled response. 

Thesis Supervisor: H. Philip Whitaker 

Title: Associate Professor of 

Aeronautics and Astironautics 
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OBJECT 



The object of this paper is to analyze a surface-to-air missile 
using equations which describe its aerodynamic and control system 
responses; to select values for control system parameters which will 
give specified response; and to investigate the possibility of compen^ 
sation which will eliminate coupled response effects between sub- 
systems. 
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CHAPTER 1 



INTRODUCTION 



1. 1 Physical Description of the Missile System 

The missile studied is a surface launched antiaircraft weapon. 

It is a homing missile using proportional navigation. Thus the missile 
borne receiving antenna tracks the target and the guidance system 
continually measures the angular rate of change of the line of sight to 
the target. The guidance computer generates command signals for 
yaw and pitch autopilot systems. The autopilot systems react by 
changing the missile heading at a rate proportional to the rate of change 
of the line of sight. Because of cruciform missile shape, the yaw and 
pitch autopilots have the same form. Polarization of the receiving 
antenna and the insertion of additional autopilot commands in earth 
coordinates require roll stabilization of the missile. 

The general control problems of the autopilot, in response to 
commands generated by the guidance computer, may be stated as: 

1. An output acceleration in yaw equal to the command 
acceleration in yaw with as fast a response as possible. 

2. Same response as above, but in pitch. 

3. Maintenance of roll stabilization (zero roll angle) in the 
presence of induced roll torques. 

1. 2 System Equations 

The coordinate system to be used is shown in Figure 1.1. 

The axes are centered at c.g. and rotate with the missile airframe. 
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Fig, 1. 1 Missile Axis Orientation 



Consider i, j, k to be unit vectors along the x,y,z axes re- 



spectively. Define the following: 

F 5 applied force s F^ + Fyl' + F^ 1c (1.1) 

V z velocity of the missile = o'! + v^ + vw1c.- (1. 2) 

uT = angular velocity of the coordinate system 

= pt + qt + ric (1. 3) 

m = mass of missile (1. 4) 

M = applied moment = M ^ + M ^ + M 1c (1. 5) 

X y z 

H = angular momentum of missile 
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Consider the missile to be a rigid body, and the axes chosen to 
be principle axes. The force and moment component equations of 
motion are: 



F 

X 


= 


m (u - rv + qw) 


(1.7) 


F 

y 


= 


m (v - pw + ru) 


(1.8) 


Fz 


= 


m (w - qu + pv) 


(1. 9) 


^x 


= 


I P 


(1.10) 


M 

y 


= 


I q - (I - I- ) rp 
y^ z X ■ 


(1. 11) 




= 


1 r - (I - I ) pq 
z X y 


(1.12) 


If we consider the gravity components of force as 


inputs to the 


, which are trimmed out in steady flight by command bias signals 
' write the applied forces and moments as: 


F 

X 


= 


T - QSC^^(M) 


(1.13) 


F 

y 


= 


QSCyO, iy , 6, M) 


(1.14) 


F 

z 


= 


QSG^(c^, ip , 6, M) 


(1.15) 


^x 


= 


QSdC^ (p, a, |3, 6, i^ , ip , M) 


(I.IG) 




= 


QSdC^ (p, a, 6, M) 


(1.17) 


M 

z 


= 


QSdC^ (^, iy , 6, M) 


(1.18) 



where: 



T 

Q 

M 

S 

C 



do 



Rocket thrust, aligned with x axis (1.19) 

1 2 

Dynamic pressure = p V (1. 20) 

Mach number (1. 21) 

Reference area, (usually body cross section) (1, 22) 

Nondimens ional drag coefficient (1. 23) 

Nondimensional yawi force coefficient (1. 24.) 
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% 


= 


Nondimensional pitch force coefficient 


(1.25) 


^1 




Nondimens ional rolling moment coefficient 


( 1 . 26) 




- 


Nondimensional pitch moment coefficient 


(1. 27) 


Cn 




Nondimensional rolling moment coefficient 


(1.28) 


d 




Reference moment arm (usually missile 
body diameter) 


(1.29) 




= 


Yaw angle of attack = arc sin 


(1.30) 


a 


- 


w 

Pitch angle of attack = arc sin 


(1. 31) 


i 

y 


= 


7 = Control surface deflection to produce 

yaw forces and moments 


(1.32) 


i 

P 


= 


Control surface deflection to produce pitch 
forces and moments 


(1. 33) 


6 


- 


Control surface deflection to produce rolling 
moment 


(1. 34) 



From Figure 1. 2 we see that 7 , 
actual control surface deflections 6 ,^ by 



i 



P 



V " ^(D ^<3) 

^ " ^(T) ■ ^(D ■ 



and 6 can be related to the 

(1. 35) 
(1.36) 

+ 60 (1.37) 
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Fig, 1. 2 Positive Deflection of Control Surfaces 



Now consider the missile to be in the following steady state 
condition: 

1. Constant speed (constant Mach) where 

V = V - n/u ^ + V ^ + wi ^ = vu + w (1. 38) 

o o o o o •'^o 

2. Missile trimmed to constant pitch angle of attack, o-^ . 

This means that 



M = QSdC = 0 

y m 

and = QSCjg = constant 

The above may be expressed by: 




(1.39) 

(1.40) 



(1. 41) 
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and 



constant 



(1.42) 



constant 



(1. 43) 



Allow motions of v, r and p as small perturbations about the 
steady state values noted above, and consider their products to be 
negligible. We find, on substituting into the equations in which v, r 
and p occur: 

ac ac ac 

% 

= m [ V (cos B) (d - V (sin a ) p + V (cos a ) r I (1. 44) 
I o ' f-/ r- Q ^ o o o 



' ac, ac, ac, ac. 



(1.45) 



9C dC dC 

“y - ^ ^ FT ' ' = - V'P% 

The partial derivatives above are stability derivatives^'^ and 
are obtained empirically for airframe and conditions specified by wind 
tunnel tests. 

Substitute the below assumptions and definitions into equations 
(1. 44) - (1. 46) to arrive at the final linearized aerodynamic performance 
equations of the missile, equations (1.64) - (1.66). 



cos 0=1, since 0 less than 15 



(1.47) 



Q = constant since V is constant and p is constant 
at steady altitude 



(1.48) 



m , I and I are constants since propellant burning 
X z 

is slow process 



p = <f) = roll rate 
q = 0 = pitch rate 



(1.49) 

(1.50) 
(1. 51) 




i 

s 



\ 

i 



r 

A 



B 



C 



E 



F 



G 



H 



L 



M 



N 



4> yaw rate 



QS 57. 3 



mV 

o 


a/3 


QS 57. 3 


ac 

_y 


mV 

o 


9 y 


QSd5r: 3i 

1 /v 

z 


dc 

n 

a/3 


QSd57*. 3. 


ac 

n 


I . 

z 


3T 


QSd 57 . 3 


ac^ 


I 

X 


W 


QSd 57. 3 


ac^ 


Ix 


a? 


QSd 57. 3 


ac^ 


I 

X 




QSd 57. 3 


ac^ 


^x 


dy 


QS 57. 3 




mV 

o 


aa 


QSd 57. 3 


ac 

n 


Iz 


aa 



Table A. 1, Appendix A, tabulates the above aerodynamic 
coefficients at three pitch angles: 

a = 0° , 12°, 24° 



(1.52) 
(1. 53) 

(1. 54) 

(1. 55) 

(1.56) 

(1.57) 
(1. 58) 

(1. 59) 

(1.60) 

(1.61) 

(1.62) 

(1.63) 



8 



i 



I 

I', 




I 



mi 



Yaw force equation: 



• • 

(A /3 + B 7 + M 6 ) = j3 - (j) sin + 4^ cos 



Yaw moment equation: 



C/3 + E 7 +N 6 = 4> - ( 




Negligible 

term 



Roll moment equation: 



G 6 + F0 + + L 7 = V 



(1.64) 



(1.65) 



( 1 . 66 ) 



The above development is based on information received from 
Convair (Pomona) and roughly parallels a more rigorous development 
of linearized aircraft performance equations outlined in Reference 1. 1, 
Sections II I through II - 8 . 

Having the aerodynamic performance equations we need now to 
consider the control systems in yaw and roll (pitch being considered 
constant in this analysis). 

The following sensors are used in the missile: 

F F 

1. Accelerometers which in effect measure — ^ and — ^ , 

m m 

or with suitable scale factors, (A0 + B 7 + M 6 ) and (Aa + B 7 + M 6 ). 

• * • 

2. Rate gyros which measure ^ , 6 . 

The accelerometers provide information for feedback 
comparison with command signals and the angular rate signals are 
used for damping. 

The sensor outputs used in this analysis are related to the 
sensed quantities as shown below: 

(A3 + B7 + M 6 )Ms) = G^(s)(A3 + B 7 + M 6 )(s) 

0 ' (s) = G 2 (s) 0 (s) 

I' (s) = 03 ( 3 ) 4 )(s) 
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(1.67) 

( 1 . 68 ) 
(1.69) 



I 




where the primed quantities are the sensor outputs and the sensors are 
quadratic of the form 



2 

(j 

G(s) = ~2 2 

S + 2^k) S + CJ 
'n n 

The control system equations for the autopilots considered in 
this analysis are the yaw control equation 

7 ^(s)= [K^(A^+ Bt+ M 6 )'(s)+ (K 2 + KgS) 4j'(s)1 (1.71) 



and the roll control equation 

6c(s) = - [ ^ + K5 + Kg s]J'(s) (1.72) 



These are integrated throqgh rate servos to control surface 
deflections 





y(s) = 


G^(s) . 
s ^0 


(1.73) 


and 




°5 • 






6(s) = 


— 6 (s) 

s c 


(1.74) 



where and Gg are of the form indicated in equation (1. 70). 

A functional diagram representation of the system is shown in 
Figure 1. 3. 

1. 3 Performance Specifications 

Consideration of the equations representing the aerodynamic 
response functions of Figure 1. 3 shows that interaction exists between 
(Aj3 + By + M6), ^ , if/ and 7 and 6 at > 0°. 

Since the system can be decoupled at = 0°, specifications for 
response at that pitch angle are to be met first. As set forth by 



10 




Fig. L 3 Missile System Functional Diagram 

Convair in stating the problem, the responses to be achieved (in 
uncoupled mode) are: 

1. Time constant of the yaw subsystem to step command 
(A0 + B 7 )^ to be less than .6 seconds. (Time constant here is defined 
as the time for the quantity to reach 6 3% of its steady state value. ) 

2. Time constant of the roll subsystem response to a step 
torque (6rp) to be less than . 5 second. (The time for <p to decay to 
less than 37% of its peak value). 

The only requirement in coupled modes (a^ = 12°, 24°) is that 
the system be stable and reasonably damped in response to step 
torque, 6rp , or step command, (A/3 + B^)^ . 
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In addition, other considerations have placed the following 
limits on and Kg; 

Kg ^ .40 (1.75) 

Kg < .021 (1.76) 

All gains (K^ Kg) are to be considered positive. 

Thus we have the system, and the performance specifications 
to be met. The remainder of this chapter will describe briefly the 
steps to be followed in the system analysis. 

1. 4. Idealized System 

Consider a system in which the servos and sensors have unity 
transfer functions: 



Gg = 1 (1.77) 

The missile with "perfect" sensors and servos will hereafter 
be called the idealized system. The idealized system will be con- 
sidered first, and then the effects of non-ideal servos and sensors 
will be discussed. It is sufficient to say that the idealized system is a 
reasonable approximation to the actual system if it can be shown that 
the servos and sensors introduce into the actual system ac^ditional 
modes of low amplitude which are far removed in frequency from the 
natural modes of the idealized system. This matter will be discussed 
in Chapter 6. For the idealized system we may consider the aero^ 
dynamic and control equations in the following form: 



(A^ +By + M6) = ^ - <f> sin + if cos 


(1.78) 


C/3 + E y +N 6 = 4* 


(1.79) 


1 •• 




GS + F0 + Hj3 i- L i y = (f> 


(1.80) 


T = -|[ K^(Aj3+ By + M6) + K 2 1 + Kgtf] 


(1. 81) 
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<1.82) 



+ Kg'^3 

1. 5 Procedure of Investigation 

The investigations of subsequent chapters will 

1. Discuss analytic methods of selecting gains to provide 
specified performance. 

2. Use analogue simulation (guided by the. inethods discussed 
above) to select a suitable set of gains. 

3. Discuss the possibility of compensating networks within 
the missile to eliminate or attenuate the effects of coupling on desired 
response. 
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CHAPTER 2 



ANALYTICAL STUDY OF THE COUPLED SYSTEM 



2. 1 General 

In Section 1 . 4 of Chapter 1 are shown the linearized differential 
equations of the idealized system. These can be put into the form of a 
"system determinant" by gathering all of the variables in each equation 
to one side, the other side being zero. If this be done, and LaPlace 
transform notation be adopted, there results the system determinant: 



0 


• 


t 

</> 


7 


6 


(^-A) 


cos 


a -sin a 

o o 


-B 


-M 


-C 


s 


0 


-E 


-N 


-H 


0 


(s-F) 


-L 


-G 


0 


0 


X 

2 


0 


1 






s 






-K^A 


-Y 


0 


(s-K^B) 


-KjM 


where 


X 


^ (Kg + KjS 


+ K4) 






Y 


= (KgS + K^) 







Expansion of this determinant yields the system characteristic 
equation. This is a sixth order polynomial in s in whose coefficients 
appear a considerably involved mixture of the controller gains. 
Appendix B gives this characteristic equation where the aerodynamic 
coefficients for = 24 ° have been entered. 

The objectives of the analytical investigation are: 

1 . To determine how to fix upon a set of gains through 
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Kg which produce a stable system and 

2. To determine a set of gains which will provide for meeting 
the specified preformance at = 0 described in Chapter 1, Section 1. 3 
and 

3. To select these gains not only to meet objectives 1. and 2, , 
but further to obtain conditions in which the disturbance in roll caused 
by insertion of yaw command and the disturbance in yaw casued by roll 
torque disturbances are acceptably low and of reasonable darrp ing. 

Inspection of the characteristic equation reveals that several 
inequalities can be obtained from the expressions for the coefficients 
of powers of s which must be satisfied to prevent occurence of roots 
of the characteristic equation in the right half of the s plane. Given 
a set of gains, the test for presence of sign change in the equation 
can be readily performed. 

It would also be possible to apply the Routh Criterion, at the 
expense of additional labor, for a given set of gains, to obtain further 
information as to the nature of the roots. 

These, and other procedures that might be derived afford little 
in the way of a guide to selection of a new set of gains to test when one 
set has resulted in an unstable condition. 

There is the possibility of forming from (2.1) a set of non- 
interacting transfer functions from which a block diagram of the system 
could be constructed. With manipulation of this diagram and the use of 
root locus methods, the roots for a given set of gains could be found. 
Unfortunately there would be at least three inner loops of such a diagram 
whose roots are gain dependent, and a very considerable amount of 
tedious and repetitious graphical work would be involved, even if each 
new attempt could be based on a set of gains choice of which was 
guided by the last trial. 

It is apparent that it would be desirable to have a more orderly 
approach in which less futile graphic or algebraic work is expended in 
arriving at the conclusion that the last guess was unsatisfactory before 
making a new, more or less intelligent, guess. 
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The object of the following section is to discuss a procedure 
which, while not avoiding trial and error and not avoiding tedium, at 
least puts to good use the experience gained by past results. 

This procedure begins with manipulations of the characteristic 
equations of the uncoupled roll and yaw subsystems. These equations 
are of low order and quite straight forward so that the effect of gain 
changes on the resultant responses are readily seen. 



2, 2 An Algebraic Approach to a Stable System Characteristic 

Equation 



In Chapter 3 is derived a block diagram interpretation of the 
system equations which results in a new form of the system determi'- 
nant where in the coupling effects and the subsystems are clearly 
depicted. 

When (3,7) o^t Chapter 3 is reformed forest^ = 0 we have 



Yaw Subsystem 
1 -F -F 



12 



-F. 



A = 



-F 



1 

-F, 



-F, 






3 " 1 

0 0 

0 0 



0 

0 



0 0 
0 0 
0 0 



l\ 

-F 



11 



10 



( 2 , 2 ) 



Roll Subsystem 



Expansion of (2, 2) will reveal a result which is the product of 
the results obtained by expanding those vsections of (2, 2) labelled Yaw 
Subsystem and Roll Subsystem as if they were themselves determiJ- 
nants. 

When the block diagram Fig, 3, 7 of Chapter 3 is considered 
at = 0 that diagram divides into two new diagrams of the roll and 
yaw subsystems. Using the methods described in Chapter 3, one can 
obtain the subsystem determinants which are the subdeterminants 
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within (2. 2). 



Then if we let 

1 -F 



A 



10 



-^11 1 



s ^ (s -F) + G(Kg + Kg e + K^) 
(s-F) 



(2.3) 



a,nd 



y 



-F 



13 



-F^ 1 



-^3 -Fi 



12 



-F. 



[ +(-K^B-A-EK 3 ) + (AK^ B-C-EK 2 + EAKg - BCK 3 - BK^A) ^ 

i(s-K^B)( 



(CK^B = EAK 2 + EAK^ - BCK 2 ) 



(s-Kj^B)(s-A) 



(2.4) 



Employing the same analog with Cramer's Rule as is used in 

Chapter 3 we can express the outputs due (A^ + B7)^ and 6^ in this 

a =0 case as 
o 



(^(s) 

6^ 



-F 



10 



A 



(2.5) 



(A3 -t B7 (s) _ 
(A3 + Bt(s) ■ 



AF 



1 


-F 

13 


1 




1 


-F 

^13 


-^^12 


-^6 


1 


0 




0 


.1 




CO 

1 


-^1 


0 




0 


-^1 


1 



( 2 . 6 ) 



When (2. 5) , (2.6) are expanded, substitutions made for the defined 
transfer functions, and simplification made, applying a step input 



(A3 + B 7 )^(s) = and 63 ,(s) = ~ . 



17 



t 




</»(s) = 



RG 



(2.7) 



(S^ + (GKg - F) S + K^) 



-K. PA ( Bs - El 



(A/3 + Bt )(s)= ^ ^ 

S [ S + (-K^B - A - EK 3 ) S + 




-K^ PB [ S (S-A)+ C] 



( 2 . 8 ) 



In the results of (2.7) and (2. 8 ) it is clear that the modes which 
exist are completely specified by the denominators and that in each 
case these are 3 in number. It is further apparent that the zeros of 
both closed loop functions are due to aerodynamic coefficients. 



characteristic equations on the assumption that the performance is 
dominated by a pair of complex roots and has prepared graphical aids 
to the rapid determination of the pole-zero configuration of the closed 
loop function which reflects a specified time-domain response to a 
step input. With these it is possible to rapidly obtain a set of gains 
for the separate subsystems to achieve a specified response. However, 
there is no certainty that the assumption of " 2 nd order dominance" in 
the subsystem will produce a set of gains which will make the coupled 
system stable. 

It is possible, though, to assume a set of gains and rapidly 

determine the residue coefficients of the modes that exist using methods 

2 2 

described by Thaler. ' By this method a check can be quickly 
obtained upon the ability of a stated set of gains to meet the specifir 
cations at = 0. A map can then be prepared of the locations of roots 
of the characteristic equations of the subsystems which do satisfy the 
a - 0 . specifications. 



Having obtained a starting set of values of the gains, these can 
be entered into the coupled system characteristic equation (Appendix B) 
and that equation factored, the roots plotted with the roots of the 
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characteristic equations of the uncoupled subsystems, and the trend 
no ted, 

A second trial for the roots can be guided by the concept that 
increased stability of the coupled system will be accompanied by 
increased "sluggishness” in the = 0 modes of the subsystems. In 
the actual carrying out of this process it is also noted that transient 
oscillating frequencies of = 0 modes tend higher from the first 
choice of subsystem root locations. By following a path on which the 
complex pair of roots of each subsystem increase in natural frequency 
but remain at constant damping ratio the real root moves toward the 
origin, and convergence to stability of the coupled system characteristic 
equation is apparently more rapid than on other paths. 

Thus there does exist one orderly, albeit tedious, process for 
converging upon a set of gains which produces stability of the coupled 
system and yet allows the subsystems to meet the = 0 specifications. 

"Bold spotting" allowed convergence to a stable set of gains in 
three tries using these concepts. 

Somewhat less tedious labor is involved if one does not observe 
the constraint of specifications too rigorously, merely working with the 
subsystem characteristic equations to produce stability in the coupled 
system characteristic equation without too much regard for specifications, 
then trimming to achieve the specifications at = 0. 

Such a procedure is, then, capable of satisfying objective (1) of 
Section 2, 1. It is quite apparent that repetitive factoring of the coupled 
system; characteristic equation as would be necessary with this (and 
most other) methods to reach the other objectives is simply beyond 
reason without automatic computing facilities. 

The conclusion is therefore drawn that this method does set forth 
a concept which is useful in guiding one to a stable set of gains for an 
analogue simulation of the problem, but is otherwise of academic interest, 
and totally devoid of practicality in achieving the other objectives of 
Section 2. 1. 

For this reason all work in "trimming up" the system was done by 
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analogue simulation as discussed in Chapter 4. 



2. 3 Steady State Analysis of the System 

In the preceding section it was concluded that setting the gains 
in the system to meet all the stated objectives is not a practical task for 
analytical study, that rather it is an undertaking more appropriate to 
the analogue computer. Considerable useful information for a com- 
puter study can be gained by analyzing the features of the variables 
considered as outputs in the steady state. 

The methods described in Chapter 3 for solving the equations 
will be used here. The uncoupled mode at = 0 will be dealt with 
first. 



1, Uncoupled Roll Subsystem 

In Section 2. 2 it was shown that at =0 the roll system charac- 
teristic equation is given by 



-F 



10 



-F 



11 



S^(s-F) + G (Kg + KgS + K^) 
(.6 - F) 



(2.3) 



It was also shown that 






(2. 5) 



Since 




ap^jlicable. For 



is analytic in the RHP, 




the final value theorem is 



Lim </)(t) = Lim [ "^<^(,s)J 
t~* cO s 0 



S [ 



R 



a (s - F) 



S^+ (GKg - F)^^+ GKgS +GK^ 
(:s - F) 



= 0 



(2.9) 
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Thus the roll system ultimately recovers from a roll step 
disturbance. 

2. Uncoupled Yaw Subsystem 

From (2. 4) the characteristic equation of the uncoupled yaw 
subsystem is 

[ s^+(-K^B-A-EKg)s^+(AK^B-EK 2 + EAK^ - BCK^ - BK^A)s 
^y " s (s ■- B)(s - A) 

(CK. B - EAK„ + EAK, - BCK.,)] 

^ s(s-K,BKs-l) ^ 

In (2, 8) for a step input (A /3+ By)^(s) = ^ 

- PAK, [ Bs - E] 

(A/3+B7HS) = ^ 

s [ s + (BC - AE)(K^ + K 2 >] 

-PBK, [ s^ - As + C] 

+ ^ (2.5) 

s [ s + (BC- AE)(K^ + K 2 )] 



Then since for a stable choice of gains this function is 
analytic in the RHP the Final Value Theorem applies and 

Lim (A3 + By )(t) = Lim [ s (A3 + By )(s)] 

00 s -■ 0 

-PK, (BC - AE) K, 

- i = -P ( 1 ) 

+ K 2 (BC-AE) + K 2 



Applying the Cramer's Rule method of solution to the yaw 
subdeterminant in (2. 2) 



-K 



1 



^ 

(A 3 + By) 



-(s) = 



-Fa (-Fg - Fg Fg) [Es(s-A)+BC] s(s-A)(s-K^B) 



( 2 . 12 ) 
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I 




or 



(Aj3 4- B7) 



(s) = 



-K^ [ Es + (BC - AE)] 

[ + ... 



(BC - AE)(K^ + K2)] 



(2.12a) 



Using (A/3 + B7) (s) = — again, and applying the Final Value 
c s 



Theorem 



Lim 4i(t) = Lim s [4'(s)l = 



t -►00 



s 0 



Ki+ K2 



(2.13) 



Thus it is apparent that both (A/3+ By)(t) and 4<(t) come to a 
constant non zero value in steady state and that 4>(t) increases without 
limit . Their senses and magnitudes are easily calculated for com- 
parison with computer results. 



3. Coupled System 

In performing the steady state analysis at > 0 it is convenient 
to let s -*• 0 in the system determinant (3.7) at the outset since that 
determinant appears in the denominator of the Cramer' s Rule solution 
for each variable of interest. In the following statements, the fact that 

where 



Lim A 
s - 0 



1 

s 



. done will be 


signified by 


[A] 

ss 


1 


*^2 


A 


0 


M 


K^B 


B 


B 


-E 


0 


-C 


0 


-N 


B 


cos 0! 

0 


1 


sin a 

0 


M 


A 


A 




A 


A 


L 

F 


0 


H 

F 


1 


G 

F 


0 


0 


0 


^4 


0 



= A_[A]^s < 2 . 14 ) 

s 



Other notation is as explained in Chapter 3. It is assumed 
throughout that a stable situation exists so that the Final Value Theorem 
is applicable. 
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R 

a. (j> due to roll disturbances 



</,(s) (2.15) 



s [ A ] 

2 *■ rr ss 

Litn <f>(t) = Lim s [ (/>(s)] = — : 

t ->’co s -* 0 — ^ [ A ] 

_3 ss 



R s [ A ] 

*■ rr ss 



[ AJ 



(2.15a) 



ss 



Evaluation of the determinants shows <j)(oo) = 0 when a roll 
disturbance is applied. 



b. <j> due to yaw command (A0 +B7) (s) = — 

c s 



Lim (/)(t) = Lim 
t -*-oo s -* 0 



s ^yr s 

<f,(s)=— 

® ^yr^ ss 

s.[ ^(s)] = _s 



(2.16) 



(2.16a) 



Again evaluation of the determinants shows <|>(«®) = 0 when a yaw 
command is applied. 



• P 

c. ijj to a yaw command (A/3+B7) (s) = — 

c s 



A ^ 
^ (s)=-I^ 



• *. .. ^ ^ t AyyJ gg 



Lim 4^ (t) = Lim s [ 4'(s)] = 
t-*oo s-0 



[ A] 



ss 



(2.17) 



(2.17a) 



t 

Evaluation of the determinants shows that i(j goes to the constant 
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iv 




(2.17b) 






-KiP 



ss K, cos a + K„ 
1 o <s 



and that i|j increases without limit as t co 



f 

d. ij) due to roll disturbance 6rr,(s) = — 



(s) = 



s ry 



(2.18) 



Lim 4^ (t) = Lim s [ 4^ (s) ] = 
t -*09 s — 0 



R s [A ] 
^ ry^ ss 



[A] 



ss 



(2.18a) 



Evaluation of the determinant shows that i|^ (®^) = 0. 

Further evaluation of the integral shows that 4^ = 0. 

This information on steady state values is of use in correlating 
the analogue computer results in Chapter 4 with the physical operation 
of the system. 
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CHAPTER 3 



COMPENSATION 



3.1 General Approach 

The general problem of preventing interaction between con- 
trolled quantities in multiple-loop systems has been discussed by 

3 1 3 2 2 1 

Tsien, * Truxal * and others. Among these Chu * has developed 

a generalized theory of multiple-loop systems including techniques for 

the determination of compensation necessary to prevent interaction. 

Chu employs methods of manipulation with the system determinant 

with the aid of block diagrams. The methods of Tsien and Chu are 

quite parallel in theory, the choice between them being a subjective 

matter. Chu' s method has been chosen for application to this 

problem since that technique is more familiar to the authors. 

A brief description of the basic elements of Cha' s approach 
to mulitple-loop control systems is included in Appendix C and no 
further detailed explanation will be given here. 

The problem of analyzing the effects of interactions in this 
system is essentially one of maintaining contact with the physical 
hardware involved. Particularly, a block diagram possesses the 
general appeal of being a reasonably tangible representation of the 
relations among the physical variables of a system, provided that 
diagram does not endeavor to lump too many quantities and inter- 
relations into a single complicated block transfer function. 

A diagram which, it is believed, will provide a reasonably 
physical view of the problem will first be derived. It employs the 
most basic set of transfer functions obtainable from the linearized 
equations. 
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Consider the set of linearized equations in the form of a 



deterr 


n inant: 










(I) 


(S-A) 


COS 


a ^ -sin a -B 

o o 


-M 


= 0 


(II) 


-C 


s 


0 -E 


-N 


= 0 


(III) 


-H 


0 


(s-F) -L 


-G 


= 0 


(IV) 


0 


0 


^ 0 


1 


= 0 








s 






(V) 


-K^A 


-Y 


0 (s-K^B) 


-K^M 


= 0 












(3.1) 


where 


(I) 


- 


yaw force equation 








(II) 


= 


yaw moment equation 








(III) 


= 


roll moment equation 








(IV) 


- 


roll control equation 








(V) 


= 


yaw control equation 






and 


X 


= 


(Kg s^ + KgS + K^) 








Y 


z 


(K2 + KgS) 







From each equation represented by a row of the determinant 
a single variable is selected to be described as a function of all the 
others. For manipulative ease, a transfer function describing each 
relation is defined: 







• 


t 

</> 


y 


6 


s - A 




- COS a 

0 


sin a 

0 


B 


M 



Equation (I) 

o - cos a 

(s-A° ■ 
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(3. 2b) 



o sin O' 

■f 



P2<s) 



' Ti?5T • 





F4<s) 



(3,2c) 
(3. 2d) 



The theorem of superposition suggests representation of /3(s) 
by the elemental block diagram. 




Fig. 3.1 



9 

4^ 






T 


6 


9 

</> 


s 




c 


E 


N 


0 



Equation (II) 

F^(s) 

f (s) = 4: Fj.(s) 



(3.3) 



(3. 3a) 



(3. 3b) 
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I 

i 



i 




A(.) = ^ = F,(S) 

i,s) . 0 




4j(s) 



Fig. 3.2 



• 


0 


7 


* 


6 


(s-F) 


H 


L 


0 


G 



Equation (III) 

f “ T5^ ^ 

!<=) = Ti^ = V"> 



(3. 3c) 
(3. 3d) 



(3.4) 

(3.4a) 

(3.4b) 

(3.4c) 
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<?(s) 



Fig. 3. 3 



6 




y 


• 


t 


1 


0 


0 


0 


C.X 

2 

s 



Equation (IV) 



-f (s) = = F^(s) 

<{> s 



# 



^ (s) 


Fn 


6(s) 




11 





Fig. 3.4 



T 




0 


# 


6 


(A/S + Bt)^ 


(s-K^B) 






Y 


K^M 





Equation (V) 



(3.5) 



(3.5a) 



(3.6) 
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i 








It will be noted that provision has been made for the insertion 
of a yaw command signal which subtracts from the sensed side 
acceleration (A0 + By + M6) as indicated by Figure 1. 3 of Chapter 1. 



X(s) = 

^<s). 

(AP+By)^ 



KiA „ , , 

1 - F,„(s) 

(s-K^B) = 



(s-K^B) = ^13^®^ 



K^M 



= F,,(s) 



(s-K^B) - ^14' 






(3.6a) 



(3.6b) 



(3.6c) 



(3.6d) 




Fig. 3. 5 



The interrelations existing among the system variables may 
then be described by a block diagram as shown in Figure 3.6. 
Provision of a node for introduction of a roll disturbance torque 6rp 
will be noted. 
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Fig. 3.6 System Block Diagram 
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The method of Chu requires translation of the block diagram of 
a system such as depicted in Figure 3.6 into a "Standardized Block 
Diagram". In this "standardized" diagram certain manipulations to 
Figure 3. 6 are made to satisfy a set of conventions which allow a 
symmetry to exist between the "standardized" diagram and the system 
determinant and which allow the writing of the system determinant 
by inspection. These are described in Appendix C. When the 
necessary adjustments are made the "standardized" block diagram 
takes on the appearance of Figure 3.7. 

In Figure 3.7 it will be observed that the manner of redrawing 
Figure 3. 6 to produce the "standardized block diagram" has resulted 
in the division of the coupled system into regions which can be referred 
to as the roll and yaw subsystems. It is readily apparent that Fg and 
Fg represent the coupling effects of yaw into roll while Fg, F^, F^ and 
Fj^^ clearly show the coupling effects of roll into yaw. Inspection of 
the transfer functions represented by Fg, F^, F^, Fg, Fg, Fj^^ immedi- 
ately show that at a = 0 the system separates into two subsystems. 

By application of the rules (Appendix C) for writing the system 
determinant from the "standardized block diagram", that determinant 
becomes: 

Yaw Subsystem Roll to Yaw Couplings 



1 


-^^13 


-1^12 1 


0 


-1^14 


CD 

1 


1 


-^5 1 


0 


-^7 


1 

GO 

1 


-^1 


1 1 
1 


-F 

^2 


-^4 


1 

CD 

1 


0 


1 

-^8 1 


1 


-^10 


0 


0 


0 1 


-F 

^11 


1 



Yaw to Roll Couplings Roll Subsystem 
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Fig. 3.7 "Standardized Block Diagram 



» 






I 









Equating A to zero and evaluating one obtains the system 
characteristic equation. Immediate proof that above represents the 
original determinant of the system can be olptained by reading each 
row of A with the appropriate substitutions. 

As a consequence of the way selected for redrawing Figure 

3.6 to produce Figure 3.7, A divides into regions illustrative of the 

subsystems and the intercouplings as noted above. The close 

3 1 

parallel between this determinant and Tsien' s ’ matrix are 

apparent. Also the clarity with which interactions are displayed in 

Figure 3. 7 and in A would seem to some degree to refute one of the 

3 3 

objections to the block diagram approach * . 

It is true that considerable effort (although simple and direct) 
has been put into construction of this form of A , but this has placed 
the designer in the position of having a graphic representation in 
Figure 3. 7 of Cramer' s Rule. Since each node output of Figure 3. 7 
represents one of the variables of interest in the problem it is merely 
necessary to follow Cramer' s Rule to assess the output at any node 
due to input at any node. 

3. 2 Objectives and Procedures 

The objectives of complete compensation will be assumed to 
be: 

1. isolation of 0 from (A^ + 

2. isolation of ip from 6rp 

The introduction of any compensation will change A , and in 
this will clearly affect all outputs. This change should not produce 
serious effect in the desired responses, (i. e. Should not produce 
an unacceptable array of closed loop poles). 

From the foregoing statements regarding the association 
between Cramer' s Rule and Figure 3.7 and A it is clear that it is 
possible to solve for 0(s) as a function of (Aj3 + B7)^ as follows: 
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(3.8) 



r FA 

^ a yr 

(A0 + Bt) A 



Where ^yj. been formed from A by replacing the IV (or <l>) 
column of A by a column having zeros in all row elements except the 

first at which unity appears (the entrance function F and the yaw 

^ ♦ 

command signal transform having been factored out of the <f> column of 
A ). 

yr' 



Likewise 



A. 



A 

(s)=^ 



(3.9) 



The desired response in yaw is given by 

A 



A. 



(A/3 + By) 



(s) = 






(3.10) 



And the normal (existing, though not particularly desired) 
reaction of ^ to 6rp is 



1 ^ 

-L(s)= ^ 



(3.11) 



Provided the missile can be prevented from yawing due to roll 
disturbance, no side acceleration will occur. 

Then the desired goals of compensation are achieved if 



1 . 

2 . 



A = A 
yr ry 



(A/3 + By) 



(s) is not altered by accomplishment 



of (1) except in the sense that A is changed, and 
that this change is not deleterious. 
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OP 



i 

1 

JV 



It is of interest to inspect both A and A 

t' yr ry 




(3.13) 



(3.14) 



As the dashed lines indicate, the expansion of ^j,y ^ 
by cofactors pivots on the circled element and suppresses the 
indicated rows and columns. Thus 






1. a compensator placed at rl, c IV will not 
A 

yi" 

2. a compensator placed at r V, cll will not 
A 

r*y 

Also of interest are A and A : 

rr yy 



appear in 



appear in 
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I 



1 




■^13 


-^12 


f 


-F 

14 


-^2 


1 






-^7 


-^3 




■^1 


1 


1 

1 


-^4 


1 

CD 


0 


-^8 


i 

(ih 


-F 

^10 


n 




n 




u 






yy 




fT\ - 














u — 


^14 


-^6 


\ 


1 


-^5 


0 


-^7 


-^3 


J 


i 


1 


-^2 


-^4 


1 

CO 


( 

1 


) 

1 


-F 

8 


1 


-^10 


0 


\ 


1 


0 


-^11 


1 



(3.15) 



(3.16) 



From these it is noted that a compensator placed at rl, cIV does 
not appear in A and a compensator placed at rV, cll does not appear 



in A 



yy 

Thus one concludes that compensation at rl, clV can possibly 

be used to make A = 0 and a compensator at rV, cll can possibly be 

used to make A^^ = 0 and that neitoer compensator will appear in 

A and A 

yy rr 

From the way in which the system determinant is constructed 
from the "standardized block diagram" it is apparent that the rl, cIV 
element of A constitutes a feedback from node IV to node I in Figure 
3. 7. Inspection of the diagram reveals that use of this path implies 

4 e 

that </) is available and that modified can be added algebraically to 
7 . Likewise rV, cll of A constitutes feed forward from node II to 
node V. ip is physically tangible. The same question of ability to 
mix the modified ijj signal with 6 arises here. If the compensating 
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networks required are realizable either exactly or approximately this 
question of mixing can be studied further. 



Installation of compensation functions - V(s) at rV, cll and W(s) 

at rl, cIV in A , where, conveniently, no element previously existed, 

results in the solutions, when A and A are equated to zero: 

yr ry 



-Ls + (LA - BH)s + (HE - LC) cos a 



(GE-LN)s + (-AGE-LMC - HNB + LAN + HME+GCB) (3.17) 




612s^ + 143s - 14.6500 
-^1. 2s + 42200) 



(3.18) 




133. 5s^ + 47. 8s - 1548 
-[ (100, 300)s + 37,570] 



(3.19) 



sino(^(s-A)[(CG-HMKps+K^(-BCG-LAN-HME+LCM+AGE] 
" (s-K^B)[(GE-LN)s+(-AGE-LMC-NHB+LAN+HME+GCB)j 



s (s-F)[ Ns+(CM+K^ME-NA-BNK^)j 
(s-K^B)[(GE-LN)s+(-AGE-LMC-NHB+LAN+HME+GCB)j 

(3..20) 



Unfortunately neither of these networks is realizable with 

3 . 

passive circuit elements exclusively. 



Upon examination of V(s) it is found that it can be divided into 




6. 71s^ _ 1. 57s 1600 

s+46 3 s+463 s+46 3 




I. 335xl0~^s^ _ . 478xl0~^s 15. 5x10 

s + , 376 s + ,376 s + , 376 



(3.18a) 



(3.19a) 
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An approximation using the last two terms of V(s) is 

-* 24 ° 

employed in the analog computer simulation with the result that 
while addition of this approximation to V(s) does not create an 
unstable system, neither does it attenuate the roll disturbance due 
to a yaw command input. Thus the indication is that the first term 
is of dominant importance. 

Thus hope of a simple approximation to V(s) which would 

attenuate the unwanted roll response is lost. By arrangement of 

the analog simulator to provide an approximation to the first term of 

V(s) above, about 30% decrease in the extremes of roll response 
24 

to yaw command is achieved, indicating that at least partial decoupling 
can be achieved. Addition of the second term as above does not 
show any appreciable improvement over the first term alone. At 
this point available facilities on the computer are exhausted. The 
actual computer techniques employed are discussed in Chapter 4. 
Figure 5. 7 shows the improvement achieved. 

Since, if isolation is to be achieved, simulation of the first 

term of V(sM is of paramount importance, it would appear that a 
■^ 24 ° , 

tandem pair of derivatives on the available i|j sensor output is 

required. In this connection it is worthy of note that the system can 

be driven unstable quite easily when more gain than specified is used 

in the double derivative term of V(s). 

The problem of introducing the modification to mix with 6 has 

3 *5 

not been examined here. Available literature ' does not indicate 
any hope of introducing the compensating signal other than to 
differentiate it and mix it with 6^ rather than 6. That this procedure 
works has been verified on the analog computer. 

It will be noted in Figures 5. 8 and 5. 9 that the desired 

« 

responses to a yaw command, (A^ + By +M6) and ijj » have not been 
adversely affected by the approximate compensation employed. 

Relative to the possibility of using an expanded version of 
W(s) similar to that employed to approximate V(s) , it will be noted 
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that W(s) is dependent on and thus the process of adjusting gains 
would involve adjustment of the compensator. This would result in a 
"tail-chase" operation and is not felt to be practicable. Computer 
facilities available do not permit examination of the effects of 
approximate simulation of W(s). 

Since the desired approximate compensation for =12° and 

= 24 are widely different, it is apparent that if the compensation 

for = 24° is installed permanently, the compensation for = 12° 

will not be proper, even in an approximate sense. Also any 

compensation left in permanently will have the effect of creating 

coupling at = 0. Then some form of logical operation is necessary. 

It is suggested that, since the = 0 and = 12 responses are not 

greatly different (See Figures 5. 1 and 5. 2), the compensator for 

a = 24° alone be installed and that this be switched in and out. 
o 

Inasmuch as maneuvers which would be likely to develop 
extreme angles of attack would probably be presaged by some violent 
excursion of the seeker tracking line relative to the airframe 
longitudinal axis, it seems conceivable that this information could be 
used as a basis for the necessary compensation switching. This 
possibility has not been examined in detail. 

In conclusion, a method has been described for discussing 
compensation in a clear manner which involves explicit relation among 
the measured quantities available. A compensation which will 
theoretically completely isolate roll output from yaw command has 
been found and has been approximately verified on the analog 
simulator. 
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CHAPTER 4 



COMPUTER SIMULATION 



4. 1 Description of the Computer Setup 

The missile system is simulated on a Reeves Electronic 
Analog Computer (REAC), Outputs are recorded on a Sanborn 4 
channel recorder. The idealized system is instrumented, and can be 
represented by equations (4.1) through (4.5) below. These are 
rearrangements of equation (1.78) through (1, 82), and provide 5 
equations in 5 unknowns with ea,ch equation written to specify the 
highest derivative present of one of the variables. The functional 



layout of the equations is shown in Figure 4. 1. * 




• # 


0(6+ 6 rp) + F(/) + H j3 + 




<t> = 


(4.1) 


0 




(4.2) 


6 = 


-(K 4 <<» + Kgc/. + Kg,<») 


0 


• • 




^ = 


(Aj3+ B 7 +M 6 ) + sin a^(f> - cos 4* 


(4.3) 


*« 

4^ = 


C^ + E 7 + N 6 


(4.4) 


II 


Kj^[(A/3+ B 7 +M6)-(Aj3 + B 7 +M 6 )^] 






* •• 






+ + Kgilj 


(4. 5) 



^Figure D. 1, Appendix D shows the actual scaled layout used. 
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4,2 



Fi.;. 4. 1 Computer Functional Layout of Idealized System 



The system of Figure 4.1 is the one on which all gain variation 
and adjustment tests are run. In addition, the compensator network of 
Figure 4. 2, representing the approximate compensator. 



-r(s) = V(s) 
4 ^ 



24 ^ 



6.71s ^ 1. 57s 

s+463 s+463 



(4.6) 



as developed in Section 3. 2 is used in testing the effectiveness of 
decoupling compensation at > 0. 

Some modifications of the layout of Appendix D are necessary 

to permit the insertion of the scaled compensator as a transfer function 

from 4^ to 6. The modifications are minor, however, and the system 

is not shown in its modified form. The approximate derivative network 

4 1 

is from Johnson, Analog Computer Techniques . * 




Fig. 4. 2 Analogue Computer Simulation of Decoupling 

Compensator 
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4. 2 Gain Adjustment Techniques 



The following gains found by the methods of Section 2. 2 to 
stabilize the system are used as a starting set in adjustment to meet 
specifications. 



^1 - 


43. 8 


(4.7) 


^2 = 


6.18 


(4.8) 


^3 = 


. 357 


(4. 9) 


K 4 = 


27. 9 


(4.10) 


^5 = 


2. 33 


(4.11) 


^6 = 


. 01985 


(4. 12) 



Roll adjustment is done by recording and observing <f>; yaw 
adjustments made on the recorded step response of (A/3+ By + M6). 

By observation of the recorded responses, cyclical gain 
adjustments in roll and yaw, and at the 3 pitch angles selected are 
made. It is found that and Kg are to be set at or near their 
maxima. Next adjustments of and Kg in roll, and Kj^ in yaw, 
are made |;o give specified stable operation of the subsystems with 
suitable damping of the coupled modes. Subsystem frequencies (and 
thus rise time) are dependent on Kj^ and Kg; damping pn K 2 and Kg . 
Thesp are the rough guides used in obtaining the final gains listed in 
Section 5.1. Section 5.2 ^hows system response to step 6^ and to 
step (A0+ By+ at the gains selected. 

4. 3 Compensator Simulation 

A scaled setup*=^i" of the compensator of Figure 4. 2 is used in 
testing the possibility of eliminating or attenuating the effects of a yaw 
command on roll output. 



^Figure D. 3, Appendix D shows correlation between analytically 
derived response of the uncoupled mode and computer response for 
a common set of gains. 

--See Figure D. 2, Appendix D. 
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Responses of the system with and without compensator are 
compared in Figure 5. 7 through 5. 9, 
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CHAPTER 5 



RESULTS 



5, 1 Gains Selected 



Final gains selected 


are: 




^1 = 


32. 0 


(5.1) 


^2 = 


6. 40 


(5.2) 


Ks = 


. 357 


(5. 3) 


K4 = 


12. 85 


(5.4) 


^5 = 


2. 33 


(5.5) 




. 019& 


(5.6) 


System Step Responses 






Figures 5, 1 through 5, 9 are responses of <j>, 

• 

ijj to step inputs; 


(A^+Bt + M6), 


6t 


1 degree 


(5.7) 


(A^+ By+MS)^ 


= 1 deg/sec 


(5.8) 
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5. 3 System Response Parameters 



Reduction of response data from Figures 5.1 through 5.6 gives 
the uncompensated system response parameters show in Table 5.1. 
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. 06 
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00 


(A^ +Bt)^ 

- 1 deg 
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12 


.48 


. 35 


. 87 


24 


. 57 


. 23 
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Table 5.1 

Response Parameters of the Uncompensated System 



Table 5. 2 is reduced data from Figures 5.7 through 5. 9 and 
gives parameters of the compensa.ted system at = 24°. 



^roll 

sec 


<f> 

^max 

deg 


vaw 

sec 


(A^,+By +M6) 

' max 

deg/ sec 


. 34 


. 4.3 


.24 


. 91 



Table 5. 2 

Response Parameters of the System as 
Compensated at = 24° 

*Roll and yaw subsystem time constants as defined in Section 1. 3. 
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CHAPTER 6 



CONCLUSIONS 



6.1 Achievement of Specifications 

The uncoupled response specifications of Section 1. 3 can be 
readily achieved. At ca ^ > 0, coupling exists which must be elimin- 
ated or damped to maintain stable, specified performance. Stable 
coupled response is obtained at the expense of more sluggish 
performance in the uncoupled mode. 

Satisfactory response in both coupled and uncoupled modes 
may be obtained, but requires computer simulation for practical 
investigation and setting of system;parameters. 

6. 2 Elimination of Coupled Response by Compensation 

Decoupling by compensation is theoretically possible. However, 
the compensating devices cannot be achieved with purely passive 
networks. Also, some physical changes within the control systems 
are necessary to provide points at which compensating signals may be 
introduced and extracted. Third, since compensation is a function of 
angle of attack, there should be some manner of varying the compensator 
characteristics with angle of attack. 

Since some succfess has been achieved with approximations of 
theoretical compensators, a more intensive investigation of decoupling 
compensators might evolve a compensator which is both effective and 
physically practical. 

6. 3 Effect of Non-Ideal Servos and Sensors 

Since the bandwidths of sensing devices and control surface 
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servos commonly employed in missiles of this type are of the order of 
40 to 50 cps it is clear that the factors of (B. 25) are quite considerably 
removed from the poles (s-A) and (s-F) and from the origin. 

The results of factoring the ideal system characteristic equation 
indicate that its roots are not fqr removed from the poles thereof for the 
gains selected by computer study in Chapter 4. Then by inference, such 
additional roots as are introduced in the non-ideal system by the 
presence of control surface servos and sensing devices wpuld be located 
relatively close to the poles introduced by those devices, and the 
residues 'associated with these roots would be small. 

Therefore, mathematical representation by the idealized system 
is considered correct and no gross performance characteristics are 
thereby over looked. Additionally, as noted elsewhere, available 
analogue computer facilities do not permit simulation of the servos and 
sensing devices. Conversation with Professor W, L. Markey of M. I. T. 
relative to the relation of additional modes introduced by non-ideal 
instruments, to the basic modes of such systems, (especially when these 
instruments have wide bandwidths compared to the natural modes of 
the system) would seem to comfirm that this simplifying action is in 
keeping with common engineering practice. 
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APPENDIX A 



AERODYNAMIC COEFFICIENTS 



Values of the aerodynamic coefficients defined by equations 
(I. 53) through (1.62) and evaluated at the selected values of 
considered in this thesis are shown below in Table A. I. The values 
were obtained by Convair (Pomona) engineering personnel. Dimen- 
sions are such as to make all angles considered to be measured in 
degrees. 
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TABLE A, 1 

AERODYNAMIC COEFFICIENTS 
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APPENDIX B 



SYSTEM CHARACTERISTIC EQUATION 



B.l Ideal System Characteristic Equation 

Expansion of the determinant (3.7) gives the ideal system 
characteristic equation in symbolic form 



Upon substitution for the defined transfer functions and obtaining a common denominator 
for all terms except the first (unity) there results: 
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Upon substituting the aerodynamic coefficients evaluated at 

a = 24*^, the C. E. becomes: 
o 

A = - s^[516 s^ + (50 + 65, 400 K^)s + 65, 400 K2] 

- s^(s+1. 237)(-s^+(- 69K3 + . 0787 - . 364) s^ - (103.5 K 2 ) 

+ (Kg s^ + Kg s + K^)[-120, 000 s + 37, 200 K^] 

+ s^ (Kg + Kg s + K^)[1590 s + (-92. 8 K^ + 511. 3)] 

+ (K„ s^ + K. s + K.)(K„ s + K„)(81, 200 s + 40,650) 

B. 2 Non-Ideal system^ characteristic equation 

To obtain the characteristic equation of the coupled system 
when non-ideal sensing devices and non-ideal control surface servos 
are employed it is useful to note that these effect only the "roll control" 
and "yaw control" equations of (3.1). 

The "roll control" equation is modified to 

6(s) =-Gg(KgS^+KgS+K^)G2 (—3)^(8) (B. 4) 

s 

The "yaw control" equation is modified to 

st(s)= G4(K^(A|3 + B7+M6)G^ + G 3(K3S + K2)>I- ] (B. 5) 

When these are case in the form employed in (3.1) the 
system determinant becomes 
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(-K^G^GjA) (-G3G4Y) 0 (s-K^G^GjB) (-KG^G^M) 

(B.6) 



where G^ through G^ are as defined in Chapter 1 Section 1. 2. This 
constitutes a modification of the transfer functions defined in Chapter 
3 based on these two equations. The following modified transfer 
functions are therefore defined to replace the corresponding transfer 
functions of Chapter 3. 
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(B.8) 



^(s) 

4^ 






(s-K^G^G^B) 



13a 




K1G4G1M 

(s-K^G^G^B) 



14a 



(B. 9) 



(B.IO) 



I 

(A/3+ By)^ 



is) 



(s-K^G^G^B) 



F 

.a.a 



(B. 11) 
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With these modifications the same block diagrams. Figures 3.6 and 3,7, 
can be used and the expanded form of (3. 7) , see (B. l),be readily modified. 
For this purpose it is convenient to regroup the terms of (B. 1) as follows: 






+ F „ f-F F - F - F F 1 
12a^ 16 3 ^2 9^ 



^ Pua''l3al®'6Pl0-P5P4‘'9-*'7^'8'^3-V7 +F8F4^6 +F10F5F3I 



^ ^na‘ -'^10-®'8*'4-®'l®'5®'l0l 






+ F F.f-FF - F - FFF +FFF1 
^^lla 14a^ ^8 3 ^9 ^r8^6^^r9 5^ 



- F F - Fr^F - F F 
15 ^2 8 7 8 



(B.12) 



Upon substitution for all defined transfer functions except F^^^^ 
through Fj^^ there results 
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Since the control surface servos and the sensors are presumed* 
to be second order devices, it is possible to assess the effefits of the 
modified transfer functions through as follows: 



Define 
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(B.14) 



(B.15) 



(B.16) 



(B.17) 



(B. 18) 



Then the modified transfer functions become 



I' _ _ _ d n X 

(s + cs + d) (s + ms + n) s 



(B.19) 



K^A 



T 



12a 



(s + gs + h) (s + as + b) 



s - K^B ^ 



(s A gs + h) (s + as + b) 



(B. 20) 



K^Ahb 



2 2 

[s (s + gs + h)(s + as + b) - K^^Bhb] 



(B. 21) 



*Telephone conversation with Mr. R. L. Geisberg,CVAC, Pomona. 
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f h 

= 

(s^ + es + f) (s^ + gs + h) 
T~~ . 2 



s(s + gs + h) (s + as + b) - K^Bhb 

2 2 
(s + gs + h)(s + as + b) 



f h (s^ + as + ^) 

2 2 2 
[ s (s + gs + h)(s + as + b) - Kj^Bhja] (s + eg + f) 

(B.22) 



K M —2~^ . ”2 ^ 

^ = (s + gs + h) (s + as + b) 

^14a ^ 2 

[ s (s + gs + h)(s + as + ^) - K^Bhb] 

2 2 
(s + gs + h)(s + as + b) 



K,Mhb 

2 2 24 ) 

[s (s + gs + h)(s + as + Jj) - Kj^Bhb] 



From these it is clear that the common denominator of equation 
(B.13) will contain the product. 

s^(s^+cs+d)(s^+ms+n)[s(s^+gs+h)(s^+as+b)-Kj^Bbt>] (s^+es+f)(s-A)(s-F) 

(B.24) 

3 

Of these poles pf A , s (s-A)(s-F) are as in the ideal system. The 
other poles are due to the factors of 

(s^+cs+d)(s^+ms+n)[s(s^+gs+h)(s^+as+ta) - Kj^Bhb] (s^+es+f)^ (B. 25) 
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APPENDIX C 



CHU' S METHOD 



It is the purpose of this appendix to briefly describe the pro- 
cedures by which the "system determinant" employed in Chapters 2 
and 3 is derived from the "standard block diagram", A complete 
development of these standardized manipulations as well as a con- 
siderable body of analysis and design theory for multiple loop 

' (C 1) (C 2) 

systems can be found in the works of Chu ' and Thaler ' , 

It is clear from the manipulation performed in the early part 
of Chapter 3 that the principle of block diagram representation of a 
set of linear integro -differential equations is simply one of expressing 
in LaPlace notation the relations in each equation as a summation 
process wherein several quantities are employed to form a single 
output quantity from a node. This process is repeated for each 
equation, and finally, interconnections between identical variables are 
made to complete the block diagram. This procedure is in common 
use. 

From the intermediate block diagram thus derived, the system 

characteristic equation can be gotten by any one of a number of 

(C 3) 

methods including that due to Mason ' as noted in Chapters 2 and 3, 

(C 1) 

however, the method developed by Chu * provides additional facility 
in performing analysis for each of the several system variables and 
in design of compensation. 
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Figure C. 1 is the block diagram of an "n node system" (n equations in n 
variables). The diagram permits of consideration of the introduction of 
any number of inputs (forcing functions) to any combination of nodes. Any 
block diagram derived by the common process stated above can be recast 
into the "standardized" form by simply listing the nodes from .left to 
right in any sequence and performing the simple diagram manipulations 
necessary to accommodate the following conventions: 

a. Each nodal point may receive any number of signals, but 
delivers only one output. 

b. Signals arriving at a node from a node to the left (feed 
forward) are arbitrarily considered to enter a node 
additively. (Sign adjustment in the feed forward transfer 
function may be necessary to accommodate this convention. ) 

c. Signals arriving at a node from a point to the right (feedback) 
are arbitrarily considered to enter a node subtractively. 

(Sign adjustment in the feedback transfer function may be 
necessary to accommodate this convention. ) 

d. Signals to be fed forward are derived at a connection 
immediately to the right of a node. 

e. Signals to be fed back are derived at a connection immedi- 
ately to the left of a node, (This and (d) frequently require 
moving a "pickoff point" where a signal actually exists in 
the basic block diagram to the right or left past a transfer 
function block. ) 

Referring to Figure C.l and considering the equilibrium equality 
at node a we have, for the sign conventions noted above: 



a = - (G . G . )a - (G. G„ ) b - -(G G„ )n 

ab Aa be ' • no Na 

+ G, I + G„ II + G iA 

la 2a na>^» 



(C.l) 
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This equation can be rewritten as; 



m 



-n 



(1 + G ^ ^TVT ~ ^ (G )m 

ab Aa be Ba no Na ma' 



m=I 



(C.la) 



At node b we have the equality; 



b = (G , )a - (G. G„. )b - (G ,G^, )c - - (G G_ )n 

ab be Bb cd Cb no Nb 



+ + G2j^II + 



nb 



(C.2) 



which can be rearranged to give: 



- “ + °bc%b )» + <°cd°Cb>^ ■ + <Gno<5Nb)" 

= 1 <°mb'“ 
m=I 



Extension of this process to several nodes using care to maintain the 
sequence established in (C.la) and (C. 2a) it can be seen that the 
sequence of nodal equations of equilibrium can be expressed in a 
determinant: 



A = 





Gk G„ 
be Ba 


^cd^Ca ' • • : 


^no^Na 




-°ab 


1+G. G„, 
be Bb 


^cd^Cb 


^no^Nb 




-G 

ac 


be 


l+^cd^Cc • • • 


G G,. 
no Nc 


(C.3) 


-G 

an 


-G, 

bn 


-^cn 


1+G G^, 

no Nn 





Inspection of (C. 3) leads to the conclusion that, given the block diagram 
in standard form it is possible to write down the elements of the system 
determinant by inspection of the diagram, using only so much of the n 
node determinant pattern as applies to the system at hand. The main 
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diagonal divides feedback and feed forward regions and terms of the 

main diagonal are the characteristic equations of the interior loops at 

the node indicated by the row and column of the main diagonal element. 

Expansion of the determinant in any given case produces the system 

characteristic equation without the often confusing mental bookkeeping 

(C 3) 

necessary with some other methods, including that of Mason ‘ . 

By suitable choice of the method of selecting the sequence of 
nodes for drawing the "standard block diagram" it is possible (as in 
Chapter 3) to graphically illustrate in the determinant the couplings 
that exist within the original equations. Solution of the equations for 
desired input - output relations follows Cramer '-s Rule, the column 
into which the forcing functions are substituted, representing the node 
in the standard block diagram at which the desired variable exists. 

The system being linear, superposition is applicable. 
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APPENDIX D 



SCALED COMPUTER SETUP AND CORRELATION RESULTS 



D. I. Scaled Layout 

The REAC layout of the simulated system is shown in Figure 
D. 1. Note that a time scale factor of 10:1 is used (i. e. , response is 
slowed so that the computer takes 10 seconds to respond as the real 
system does in 1 second). 

D. 2. Scaled Compensator 

The compensator V(s) of Chapter 3 is approximated as shown 
in Figure D. 2. It is limited in both the inability to produce a perfect 
derivative and by some approximation of time constants to conform to 
the RC values available. 




Fig. D. 2 Scaled Compensator Layout 
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Fig. D, I Scaled REAC Layout of Idealized System 



D. 3 Response Correlation 



Figures D. 3 and D, 4 are correlation check runs between 
computer and analytic responses of the uncoupled roll and yaw sub- 
systems. 

Gains used are not the final selected values of equations (5.1) 
through (5. 6) but are listed below: 



^1 = 


43. 8 


(D.l) 




6. 18 


(D.2) 


^3 = 


, 357 


(D. 3) 




27. 9 


(D. 4) 


^5 = 


2. 33 


(D. 5) 


^6 = 


. 01985 


(D.6) 



When these values, analytic responses below are obtained: 

= R[.464e“^^*'^®^ + .465e sin(56, 28t-85. 5°)] 

(D.7) 



and 



4<(t) = -P[. 840 + . 793e“^' +. 907e”’^’ ^^^sin(16. 63t-106. 4°)] 



where 



P 5 



(A/3+ By ) I = 1 ^ 
^ c| sec 



(D. 8) 
(D. 9) 



and 

R= |6^| =ldeg (D.IO) 

In Figures D. 3 and D. 4, the continuous lines are computer 
response, circled points are check points from equations (D. 7 and (D. 8). 
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